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Abstract. This paper describes the construction of realistic 
extragalactic Monte-Carlo catalogues, aimed at comparing the 
behaviour of cosmological tests versus input parameters. QSO 
catalogues are built on a luminosity function derived from data 
through a certain choice of any cosmological and evolution- 
ary model, and through suitable computation of individual 
maximum volumes in complete (but magnitude- and redshift- 
limited) samples. Details of this computation are examined for 
the case of strong luminosity evolution. The values of the evolu- 
tion parameter are derived for various cosmologies, correspond- 
ing to < V/Vmax >= 1/2 in the sample of 400 Ultra- Violet 
Excess (UVX) QSOs (?). The various luminosity functions are 
compared, both for the whole sample and in redshift bins. A 
characteristic evolution time is defined and computed, depend- 
ing strongly on the cosmology, but practically constant when 
expressed in terms of the age of the universe. Algorithms are 
given for producing catalogues based on the null hypothesis 
that the objects are uniformly distributed in volume but suffer 
pure luminosity evolution. The only input based on real data 
is the luminosity function, but this requires neither a redshift 
nor an apparent magnitude histogram. Simulated redshift and 
absolute magnitude histogrames are compared to real data. 
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1. Introduction 

Various judicious tests have been developed to discriminate 
between different cosmological models, all requiring large and 
complete samples of extragalactic objects (quasars or galax- 
ies). Unfortunately, the available samples are magnitude lim- 
ited, selection effects induce poorly known systematic biases in 
magnitude and redshift, and extragalactic populations undergo 
strong evolution over a Hubble time. This paper describes the 
construction of extragalactic synthetic catalogues with various 
parent cosmologies, evolution laws and simulated observational 
biases, with the goal of applying them to the cosmological test 
described in Paper I (?). A preliminary study of the luminosity 
function of the real population and of its evolution in various 
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cosmologies is necessary for such a construction; so it is also a 
by product of the present paper. 

The two main extragalactic populations, quasars and 
galaxies, may be simulated: we focus on quasars in the fol- 
lowing since they are more distant and a priori more suited 
to cosmological tests. A quasar being by definition stellar in 
appearance, the pair of apparent magnitude-redshift (m, z) is 
sufficient to characterize it with respect to current cosmologi- 
cal tests. The result of the simulation is thus a list of apparent 
magnitudes and redshifts, not an image. 

The observed population suffers strong evolution, as 
evinced by either the standard < V/Vmax > analyse or di- 
rect derivation of the luminosity function in redshift bins. Pure 
Luminosity Evolution (PLE) seems to be the best model (?), 
but Pure Density Evolution (PDE) or mixed models are not 
rejected; and the question remains open for redshifts greater 
than 2.2 (?). Paper I describes a cosmological test which ac- 
counts for such evolution. 

Quasars also suffer various biases in their redshift distri- 
bution and a strong variability in luminosity. Our initial aim 
was to construct synthetic quasar catalogues capable of quan- 
tifying the behaviour of the test as a function of the various 
parameters: 

— sensitivity to the cosmological parameters; 

— photometric uncertainties and strong quasar variability; 

— redshift uncertainties and biasing in certain redshift ranges; 

— catalogue size and redshift limitations, leading to a pre- 
ferred observational strategy for the construction of quasar 
samples; 

— this list is certainly not exhaustive. 

Last, but not least, most cosmological tests using QSOs rely on 
some basic assumption concerning their evolution. As a conse- 
quence, their result in principle applies only within this frame- 
work. With Monte-Carlo catalogues, it is possible to test the 
uncertainty induced by any false hypothesis on evolution by 
making cross-tests between different assumed evolutions: this 
may concern the nature of the evolution, either Density or Lu- 
minosity Evolution (DE/LE), as well as its functional form, 
PoWer law or EXPonential (PWLE/LEXP). The present pa- 
per addresses only the case of PLE. PDE will be the subject 
of a forthcoming paper. 

The paper is organized as follows: in Section 2 the basic 
method is exposed. After a brief review of the cosmological 



2 



frame, the computation of volumes and maximum volumes is 
detailed in Section 3. In Section 4 the sample of Boyle is used 
to derive the necessary LF. In Section 5 the evolutionary pa- 
rameters and characteristic times found in various cosmologi- 
cal models are compared. In Section 6 a general algorithm to 
construct artificial QSO samples is derived from our working 
hypotheses. The redshift and magnitude histogrames of syn- 
thetic catalogues are examined. 

2. Basic Method 

Previous numerical simulations have been performed in a sim- 
ilar way, essentially aimed at producing catalogues of galax- 
ies. Chokshi et al. (?) simulate galaxies in clusters uniformly 
distributed in space. In each cluster, about 1,000 galaxies are 
drawn, all at the same redshift. The magnitudes are sorted ac- 
cording to the Schechter Luminosity Function (LF), brighter 
than a certain cut-off which of course depends on the cluster 
redshift. Yee (?) first sorts the apparent magnitude according 
to the observed distribution, the redshift is then computed in 
bins of apparent magnitude according to a certain Luminosity 
Function (LF) based on a given, fixed, cosmological and evo- 
lutionary model. Our method is slightly different: first, a func- 
tional form for the evolution is chosen, and the cosmological 
parameters are varied on a certain grid. There is then a unique 
evolution parameter ensuring the condition < V/Va >= 1/2 
for a given, actual, sample, Va being the available comoving 
volume according to Avni et al (?). The zero-redshift LF $ is 
fixed by such a procedure, and fitted according to some pre- 
vious analytic choice for the LF. Once the evolution of the 
LF is known, V'a(A/o) may be computed in terms of the zero- 
redshift absolute magnitude (Mo). The resulting probability 
distribution function in the sample of absolute magnitude is 
then $(Mo) x Va(Mo). It is easy to sort the absolute magni- 
tude of object i according to this PDF. A second parameter, 
either the apparent magnitude or the redshift, is needed to 
define a QSO. Given the absolute magnitude Mo and the cos- 
mological and evolutionary model, Va(Mo) may be computed. 
A second random drawing according to a uniform distribution 
gives a number x £ [0, 1], and the redshift of object i may be 
computed from the condition: V = x x Va(Mo). (m,z), the 
apparent magnitude-redshift pair, is then derived from these 
two variables. For each cosmological model, the procedure is 
repeated until the desired catalogue size is reached. 

The null hypotheses are: that PLE applies, that the lu- 
minosities are distributed according to the observed LF and 
that, according to the Cosmological Principle, the QSOs are 
uniformly distributed in volume on Gpc scales, so that the 
variable V/Va is uniformly spread over [0, 1]. From the LF we 
are able to draw the two variables Mo, the zero-redshift abso- 
lute magnitude, and x; and from the assumption of uniformity 
we draw the V/Vmax ratio. Within this framework, particular 
attention must be paid to the computation of individual V/Va- 

3. Theoretical background 
3.1. Cosmology 

The universe is described by a Friedman-Robertson- Walker 
model which is defined by three independent fundamental pa- 
rameters (calculated at z = 0): the matter density flrn, which 
may include a non-baryonic component, the cosmological con- 
stant JIa, and the Hubble parameter Hq. The value of Hq is 



very important for the creation of synthetic catalogues because 
it scales the distance and the absolute luminosity of the objects. 
Nevertheless, Ho has no influence on individual V/V^s, which 
are scale free. The other known cosmological parameters are 
related to these three parameters. 

No arbitrary choice of a preferred cosmological model is 
made. General formulae for the distance and the volume are 
given in CaroU et al. (?). We deflne the conformal coordinate 
X{z) as: 

Xiz) = 21^^!'/' X 

X {{i + z'y{i + nmz')- z'{2 + z')nAV/^ ' 

The volume out to redshift z is written: 

^(^) = (i)'2?v^'^^""(^(^))-^(^)i' 

where ^Ik ~ l~^m—^A is the curvature, sinn is sin if < 
and sink if > 0. This equation is very useful for numerical 
computation because it does not require inverse trigonometric 
functions. Contrary to what happens in a closed universe, with 
the expressions given, e.g., in CaroU et al. (?), Eq.(^) is single 
valued and there is no discontinuity at the value x = as 
shown in Mathez et al. (?). 

3.2. Evolution 

The quasar population resides in an extremely young universe, 
and it has evolved up until the present epoch. It is likely that 
individual QSOs do not live for a long period of time, so we 
speak rather about the evolution in a statistical sense. In any 
case, assuming a PLE model for quasars amounts to assuming 
that the fraction of active quasars is constant with epoch, their 
luminosity decreasing with increasing time, and that they are 
in sufficient number so that we may apply this average evolu- 
tionary law to individual quasars in order to derive maximum 
redshifts. The absolute magnitude at redshift z is given by: 

M(z) = Afo -2.5Zogio(e(z)), (2) 

where Mo is the absolute magnitude at z = 0, and e{z) is the 
evolutionary law: 

e(z) = (l + z)'=- (3) 

for a Power Law Luminosity Evolution (PWLE), and 

e(z) = exp {kLt{z)/Ho) = exp {t{z)/T) (4) 

for an Exponential Luminosity Evolution (LEXP); fei is the 
evolutionary parameter, t[z) the look-back time and r the char- 
acteristic evolutionary time. Strictly speaking, as quoted by 
Bigot and Triay (?), this evolution law has only a statistical 
sense and does not apply to individual QSOs. However, to com- 
pute each V/Va value, we fictitiously move each QSO toward 
high redshifts, and it is necessary to apply this evolution law 
to each QSO. 

Assuming a power law spectrum = with a spec- 
tral index a = 0.5, the Mattig relation relates the observables 
(m, z) of a quasar to Mo, its absolute magnitude at z = 0: 
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Fig. 1. Four possible apparent magnitude — redshift plots: monotonicaly increasing, one maximum, both one maximum and one mini- 
mum, and monotonicaly decreasing. The evolutionary tracks are drawn according to the cosmology and exponential luminosity evolution 
as described by the top labels. The m = m{z) curves (solid lines) correspond to bins of constant width in absolute magnitude. The hori- 
zontal dashed line represents a constant limiting magnitude B = 21; the dashed m{z) line corresponds to the faint separator 52 between 
magnitude- limited and redshift- limited quasars (see text). Depending on the parameters (cosmology and evolution), the same quasar may 
be either magnitude- or redshift-limited. In models where m,{z) shows a maximum (for example in Fig. lb), for quasars lying between S 
and B = 21 there exists a redshift slice where m(z) is above B = 21. The corresponding volume must be subtracted from the available 
volume V2 — Vi. Quasars of the complete Boyle et al. (1990) sample are superimposed. 
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w-Mo {z) = Mo - 5 + 2.5logi 



e{z) 



(5) 



here K(z) = (1 + z)'""^-* is the K-corroction. 

The evolutionaxy tracks m{z) on the Hubble diagram are 
shown in Fig. 1 for a given set of quasars in four different 
cosmologies and evolutionary histories. The curves appear to 
be quite different. In particular the quasars in Figl.b,c,d be- 
come brighter and brighter as the redshift increases, because 
the evolution dominates the cosmological dimming at large z. 

3.3. Distribution of the QSOs in the Hubble plane 

Available complete QSO samples usually have a well defined 
limiting magnitude, towards faint magnitudes. Due to satura- 
tion effects (in the photometry or in the spectroscopy), a bright 
limiting magnitude may also apply. 

All UVX QSOs are limited to redshifts z < 2.2, precisely 
because of their selection as objects with a UV excess, a stan- 
dard way to detect quasars in optical astronomy. As the red- 
shift decreases (below z = 0.3), there is another limitation be- 
cause quasars become spatially resolved and can no longer be 
classified as stellar objects. These redshift limitations prevent 
the use of maximum (respectively minimum) redshifts outside 
the limiting redshifts of the sample. Assume that we have a 
QSO sample in the magnitude and redshift ranges [mi, 7712] 
and [21,2:2], respectively. For UVX samples we would have 
zi =0.3 and 22 =2.2. 

We define 3>(M, z) as the LF at redshift z and at absolute 
magnitude M. However, for the case of PLE it can be written 
$(M(2)); then all luminosity functions at different redshifts 
differ by an additive constant along the absolute magnitude 
axis. This allows us to work with the LF defined at a fixed 
epoch {z = 0) where $(M(2 = 0)) = 4>{Mo), whatever the 
evolutionary law (?). The number of objects in the range of 
absolute luminosity [Mo, Mo + dMo] and redshift [z,z + dz] 
writes: 



d N = (j){Mo)dMoH{m2-m)H{m-mi) x 

H(z2 — z) H(z — zi) dz , 
dz 

with the following notations: 

— dV is the comoving volume element of the shell of redshift 
z to z + dz in the solid angle lo of the survey; 

— 4>{Mo) is defined above; 

— The four Heaviside functions are the selection function, en- 
suring that apparent magnitudes and redshifts remain in- 
side 'box S', or [mi,m2] ® [zi,Z2], the limitations due to 
the selection criteria. 

We now define "D as the domain in the Hubble plane where 

all quasars are redshift limited. 

Let ztop (respectively ztot) be the redshift in the range 
[21 , 22] at which the evolutionary track passes through the max- 
imum (minimum) of the m{z) curve. Ztop (ztot) is the same for 
all quasars, since their m{z) curves are parallel, differing by an 
additive constant equal to the absolute magnitude. Consider 
a given absolute magnitude Mo, and let m-Mctop (wiMo.bot) be 
the corresponding faintest (brightest) apparent magnitude in 
box B: 

mMo,top = sup^g[^j^^2] [mMo(^;)] =mMo{ztop); 



mMo,hot = illfz€[zi,Z2] ["^A^oC^)] = 'mMoiZhot)- 



(6) 



nfiMoiz) is the apparent magnitude corresponding to the ab- 
solute magnitude Mo at redshift z. ztop is equal to 22 in the 
case of little or no luminosity evolution. Consider the m{z) 
curve corresponding to rritop = m2: it is unique; let us call it 
the 'faint separator' 52. There is also a 'bright separator' the 
unique curve >Si crossing {zbot,mbot = mi). 52 is drawn 011 Fig. 
1. Domain V is simply the domain located below S2 and above 
5i . Let Mq' be the two absolute magnitudes corresponding to 
the curves Si, (i = 1,2). The domain T> simply consists of that 
part of the Hubble plane brighter than M^'^ and fainter than 
Mq^ . We call 2?* the complement of "D with respect to the box 
B. We call V and V* quasars the quasars located in each of 
these domains. 

3.4. The V/Vmax test 

Since our method relies on V/Va values, it is necessary to re- 
view the V/Va computation, particulary in the case of strong 
evolution. Consider a quasar in the sample with (m, 2) as its 
representative point in the Hubble diagram, and define the fol- 
lowing quantities: 

— Mo: the quasar absolute magnitude at epoch 2 = 0; 

— Zmax{Mo,mi) {i = 1,2): the redshift(s) where the quasar 

has an apparent magnitude m^; 

— Vmax{Mo,mi): the corresponding comoving volume(s); 

— Va'. the available volume; 

— Vi {i = 1,2): the comoving volumes enclosed by the two 
limiting redshifts Zi. 

The various possible shapes of the evolutionary track (Fig. 1) 
are crucial for understanding the V/Va test. The quasar may 
leave the sample in one of two ways: either its magnitude leaves 
the range [mi, m2l, or its redshift leaves the range [21, 22]. The 
available volume Va does not reduce to Vmax{Mo,m2): it is the 
total volume Vz — Vi, minus ^ Voxit, the sum of the volumes 
of all slices of the universe where the evolutionary track is 
outside box B, so that it prevents the quasar from remaining 
in the sample. 

Assuming that volumes are randomly distributed over the 

range [Vi,V2], the variable uniformly distributed over [0,1] is 
the ratio x defined below. The number of quasars in the redshift 
range [21, 22] is: 



N 



(Mo)2 

(Mo)l 
(Mo)2 



dMo*(Mo) 



!"/(z2>Zmax(Afo,m2)) 



dz 



sup(zi ,Zmax(Mo ,mi) 



dV 
dz 



J(Mo)i 



dMo*(Mo)Va(Mo), 



where: 



Va{Mo) = inf{V2,Vmax{Mo,m2)) - 

SUp{Vl,Vmax{Mo,mi)) - y^^Vout- 



(7) 



(8) 



Note that for V quasars inf{z2,Zmax) = Z2 and 

sup{z-i_, Zmin) = Zi Simultaneously. 

It is straightforward to find the new variable x: 



■ = {V-Vs 



sup 



Vout)/{Vinf - Vsup -Y,^' 



out ) ■ 



(9) 
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Indeed, one shows that the mean of the ratio x over the sample 
is 1/2: 



<x> 



N 



{Mo)2 



dMo$(Mo) X 



(Mo)i 



= 1/2, V$,V(Mo)i,V(A.fo)2; 
moreover, since dV = T4(Mo)da;, we have: 

dN = dMo${Mo)Va{Mo) dx 

J{Mo)i 

= N dx, 

which characterizes a PDF of x which is uniform over [0, 1]. For 
a redshift-limited sample over [zi, Z2], taking Vsup = overes- 
timates V/Va and taking Vinf = Va underestimates V/Va- In 
either case, < x > cannot be equal to 1/2 even tn the correct 
cosmology and with the correct evolutionary law. 

Varying the shape of the evolutionary track changes the ex- 
pression of the variable x. There may exist zero, one or several 
maxima of m{z), or the quasar may become fainter than m2 at 
low redshift; in fact there may be several roots of the equation 
m{z) = 7712. Fig. 1 shows various curves m{z) for various cos- 
mologies and evolutionary histories, in order to illustrate the 
existence of redshift ranges, lying inside the range [zi, Z2], but 
for which m{z) is not in the range [mi, 7712]. Such ranges must 
be excluded from the volume computation. 

4. Luminosity Functions in real samples 

The LF is not a direct observable - on the contrary it depends 
on the choice of the cosmological model and of the evolutionary 
law. For our purpose, it is necessary to be realistic in reproduc- 
ing the data; in particular, the artificial catalogues must have 
the same LF as the real sample whatever the choice of their 
parent cosmological model. 

As discussed in Kassiola and Mathez, (?), given the cos- 
mology and the evolution, there are two ways to derive the 
LF from data: firstly one can construct the Global Luminosity 
Function (GLF), by shifting the absolute magnitudes of the 
whole sample to the present epoch z = 0, according to the 
evolution law, and secondly one can compute the Restricted 
Luminosity Functions (RLFs), i.e. the LFs in different redshift 
bins. A necessary condition for PLE to apply is that the LF 
in all redshift bins differ by a simple shift along the absolute 
magnitude axis. To do this, one has to compare the RLFs. As- 
suming that PLE is correct, the reason for choosing the GLF 
to construct the catalogues is to avoid the noise created by 
binning the data in redshift. In Section 4.3, we test the PLE 
hypothesis by comparing of the RLF and the GLF. 

4.1. Real data: the Boyle et al. sample 

The complete sample of UVX quasars of Boyle et al. (1990) is 
used. It contains 383 quasars, if we exclude the Narrow Emis- 
sion Lines (NL) and the redshifts outside the range [0.3, 2.2]. 



4.2. The Global Luminosity Function and its derivation 

We now consider the computation of the GLF for any redshift- 
and magnitude limited-QSO sample. The sample is binned in 
absolute magnitude Mo with bin width AA/q. Let Ni be the 
QSO number in the j"* absolute magnitude bin. The GLF in 
bin i is given by: 



Mo 



AMo , AMo 
— - — ,Mo A — 



with an error on 



(10) 



(11) 



(12) 



For the case of the Boyle sample, none of the eight fields in- 
tersect; thus, the available volume follows directly from Eq. 
(8): 




Va = 



k 



^inf{Vj^ax,V2) - sup{V^i„, Vi) - ^ T4ui 



(13) 



where Wfc is the solid angle of the fc*'' field. 

The GLF looks similar for all of the cosmologies we tried. 
Special attention has been paid to some of the cosmological 
models (e.g. ^ = 0, A = 0; n = 1, A = 0; ^ = 0, A = 0.8), by 
various authors which found that the LF is well fitted by either 
a single (?) or a double (?) power-law. We adopted the double 
power law-model in all cosmologies and for the two functional 
forms of evolution: 



CpmodeliMo) ■^QO.i(Mo-M,)(a + l) _|_ ]^Q0.4{Mo - M* ) O + l) 



(14) 



Af* is the characteristic absolute magnitude of the two power- 
law distribution, i.e. the knee of the distribution. Absolute 
magnitudes for the quasars are computed according to Eq. (^. 
The evolution parameter kL is computed by setting < x >= 
1/2 for the whole sample. Such a procedure differs from cur- 
rent work (?; ?), in which kL is fitted like the other parameters. 
The free parameters (j>* , Af*, q, and /3 are derived by a min- 
imization of the binned GLF: 



(15) 



where Ni,in is the number of absolute magnitude bins. An im- 
portant difference with other LF computations is that we have 
only four free parameters. 

Table 1 and Figs. 2 and 3 show a panel of GLFs computed 
for four different cosmologies, scaled to Ho — 50 kms~^ Mpc~^ . 
The general behaviour of the GLF is identical for the exponen- 
tial and power-law evolutions. The GLF is enlarged and QSOs 
are brighter for a A-dominated universe. This is due to the 
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Fig. 2. Global Luminosity Functions in four cosmologies, for a 
power— law luminosity evolution. The models a, b, c, d of Table 1 
progress from the larger to the thicker line. 



Fig. 3. Same as above, exponential luminosity evolution. The mod- 
els e, f,g, h of Table 1 progress from the larger to the thicker line. 



fact that, at fixed redshift, the distance in this type of model 
is greater than in other models, and therefore the absolute lu- 
minosities are larger. The effect is in the opposite sense for 
matter-dominated universes where the objects become fainter 
and the GLF narrower. The knee of the luminosity distribu- 
tion evolves according to the preceding remarks. The faintest, 
shortest and densest GLF is obtained for Q = 1 and A = 0. 
The evolution becomes weaker in A-dominated universes but 
has a more standard value (fe_L around 3-4 in PWLE models) 
if n is greater than 0.4, regardless of the value of A. 



Despite our different LP computation method, we find re- 
sults quite similar to those of Malhotra et al. (?) and Boyle et 
al. (?). 

4.3. The Restricted Luminosity Functions 

Restricted LPs (RLPs) have been computed from the Boyle 
sample in three different redshift bins ([0.3, f.O], [1.0, 1.7], 
[1.7,2.2]). We are interested in checking whether or not the 
PLE hypothesis is confirmed. If this is the case, any RLP may 
be deduced from the GLP by a simple shift along the absolute 
magnitude axis according to the evolution law (Eq. |^. Since 
we have computed the GLP in the last Section, we are able to 
compare our RLP to the GLP. This is done by a -procedure 
between the GLP and the three RLPs: 



E 



OGLF 



•t>RLF 



Here, 



and 



(16) 



refer to the RLP, computed from the 



data in the limited redshift bins, and (f>GLF is shifted to this 
redshift range. 

The result is shown in Table 1 for four cosmologies and 
the two evolutionary models. The probabilities P(> x^) 
high, and, except for one case (LEXP, flmat = A = 1), the 
test succeeds to better than the 15% confidence level. This 
result is in condradiction to those of a similar test done by 
Kassiola and Mathez (1991), which found strong discrepancies 
between the GLP and the RLP. However, they interpret these 
discrepancies as the lack of homogeneity of their composite 
catalogues. Since our results are consistent with PLE (favoured 
by Boyle et al., 1987), we use this hypothesis to construct our 
synthetic catalogues. 

5. Evolution characteristic time in various cosmologies 

In the exponential model, the characteristic evolution time r 
enters explicitely in Eq. (jll), the expression for e{z). This is not 
the case for the power law model - but a characteristic time 
may be defined as the look-back time to which all luminosities 
were higher by a factor e: 



= t {zk. = exp{l/kL) - 1) . 



(17) 



Table 2 (3) gives the evolutionary parameters kh for a grid 
of cosmological parameters in PWLE (LEXP). Table 4 (5) 
gives the corresponding PWLE (LEXP) characteristic times, 
computed according to Eqs. ( |r^ ) and (|j), respectively. These 
times are in units of the Hubble time H^^ . As is well known, 
the quasar evolution time is of the order of 1/(^^/(7. 7± 3.5) for 
exponential evolution (?; ?; ?), strongly dependent on the cos- 
mological model. Por power-law evolution, the characteristic 
time is //o"V(3.7±0.6). 

Tables 6 and 7 give the same characteristic times as Tables 
4 and 5, but in units of the age of the universe in each cosmo- 
logical model. With the exception of the 57 = 0, A = 1 model, 
all ratios in Table 6 are close to 0.30±0.02, implying that the 
characteristic evolution time is about equal to 1/(3. 2±0. 2) of 
the age of the universe, whatever the cosmological model, in 
PWLE models. In LEXP models, this ratio is l/(6.7±0.9). The 
characteristic times in Tables 4 and 5, apparently quite depen- 
dent of the cosmology, are indeed surprinsingly constant when 
expressed in terms of the age of the universe. 
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6. An algorithm for producing quasar catalogues 

6.1. The algorithm 

Fixing the evolutionary law, the cosmology (including the 
choice of Ho) and the catalogue limits (i.e. mi , m2 , zi , Z2 which 
allow the box B and domains V and T>* to be defined), each 
QSO (m, z) is fully determined by the other pair of variables 
iy/Va, Mo). This is easy to understand from Fig. 1, where the 
Mattig functions appear. All of our previous discussion shows 
that under the PLE hypothesis (which was checked on the 
data), any sample is fully determined from its limiting mag- 
nitudes and redshifts, the GLF and a uniform distribution of 
the V/Va. From the probability distribution function (PDF) of 
absolute magnitudes Ma, which is 4>{Mo) x Va{Mo), and from 
the V/Va{Mo) ratios uniformly spread over [0, 1], it is possible 
to draw Monte-Carlo catalogues. 

The algorithm is shown in Fig. 4. Choosing the cosmology 
and a functional form for the evolution, the first part of the 
work is to compute the evolutionary parameter from the con- 
straint < V/Va >= 0.5, using the Boyle sample. At the same 
time, Va{Mo) is computed and used to obtain the GLF. In the 
second part, two series of random numbers, ei and €2, axe used 
to extract the V/Va and the absolute magnitude Mo of each 
object in the synthetic catalogue from the relations: 

V — Vsup — Vput 

£1 = V/Va = — — _^ , 

rinf y sup / j y out 

J^'^dM(l>{M)Va{M) 
J^^ dM<j>{M)VaiM)' 

Obtaining the redshifts and apparent magnitudes is 
straightforward from a simple inversion procedure. Synthetic 
catalogues have been produced in the past, where both a LF 
plus a redshift histogram were jointly fitted. As already ex- 
plained, this method suffers from at least two problems: Es- 
sentially, the two distributions may be incompatible in differ- 
ent cosmologies. Our method has the advantage of avoiding 
these problems since we use only one distribution, the GLF, 
and the other distributions (redshift, magnitude) are derived 
only after the cosmology and the evolution have been fixed. The 
algorithm shown in Fig. 4 ensures the coherence of these distri- 
butions, moreover it takes the evolution into account very ac- 
curately. This is not obviously true in previous methods where 
it is not clear how the evolution intervenes. In fact, we used 
the Cosmological Principle to replace the redshift distribution 
of the standard method by demanding the uniformity of V/Va- 
Moreover, the limiting magnitude and the redshift range for 
the synthetic sample may be different from their values in the 
input sample. 

6.2. Redshift and Magnitude Histogrames 

The input (?) and output (synthetic catalogue with the same 
observational biases) redshift histogrames are compared in Fig. 
5 for models b and / of Table 1. Similarly, the histogrames of 
absolute magnitude are compared in Fig. 6. The fit is quite 
satisfactory for both models; this is true for all models we 
tried. It is possible to introduce various biases in the synthetic 
catalogues, depending on which sort of catalogue we want to 
simulate. For example, for UVX quasars it is well known that 



Fig. 4. Algorithm for the construction of synthetic catalogues 

the color selection criteria introduce a deficiency of QSOs in 
the redshift range [1.5, 1.8] (?). Moreover, QSOs arc variable 
objects, which can be simulated by randomly selecting a mag- 
nitude variation in a pre-defined distribution function. Noise 
in the magnitude and redshift measurements may be simulated 
in the same way. 

7. Conclusion 

We have reviewed in detail the V/Vmax computation in various 
cosmologies. The Mattig relation is used to understand how to 
compute the available volume Va of a QSO, even in complicated 
cases. Both the Pure Luminosity Evolution (PLE) assumption 
and the Cosmological Principle allow the variables (Mq, V/Va) 
to be used instead of (m, z) to construct synthetic QSO samples 
in a coherent way. Some of the previous catalogue construc- 
tions were shown to suffer problems. The Global Luminosity 
Function used for this purpose is derived in a slightly different 
way from previous calculations since we fixed the evolution pa- 
rameter by the condition < V/Va >= 1/2. This significantly 
changes the value of this parameter, and hence of the GLF. 

A surprising result is that the QSO characteristic evolu- 
tion times arc constant when expressed in terms of the age of 
the universe, regardless of the cosmology. The PLE hypothesis 
w£is checked on the data, which comfirms previous results (?). 
The advantage of such catalogues is that we tightly control all 
aspects of the sample: the parent cosmology, the strength and 
the functional form of the evolution, the luminosity function, 
the magnitude and redshift depths, and the effects of magni- 
tude and redshift biases, all of which infiuence the results of 
the new cosmological test introduced in Paper I. In Paper I 
we apply the test to catalogues with different redshift limita- 



Fig. 5. Rcdshift Histogrames; (1): Boyle sample; (2): simulated cat- 
alogue, translated 15 upwards, model b; (3): idem, translated 35 
upwards, model / 



Fig. 7. Absolute Magnitude Histogrames: dashed curve: Boyle sam- 
ple; full line: simulated catalogue, model / 



Table 2. Power Law Luminosity Evolution (PWLE) parame- 
ter in different cosmologies. The model f2 = 0, A = 1, which 
is marked with an asterisk, is close to a universe without a 
Big-Bang, and not very reliable. 



Fig. 6. Absolute Magnitude Histogrames: dashed curve: Boyle sam- 
ple; full line: simulated catalogue, model b 
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Table 3. Same as Table 2 for LEXP, Exponential Luminosity 
Evolution. 



tions, which leads us to propose a new observational strategy 
for the construction of future QSO samples. It is possible to 
test a catalogue with a wrong hypothesis, for example on the 
evolution, by assuming a power-law evolution and running the 
cosmological test with the exponential hypothesis. Moreover, 
our method is tractable with a Pure Density Evolution hypoth- 
esis, provided that density-weighted volume elements pdV are 
substituted for all volume elements dV. 
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Table 1. The LF (Eq. (|l^) in various redshift ranges for four arbitrary cosmologies and two evolution laws. Since we are 
interested in input parameters to construct catalogues, confidence intervals are not needed. Only the best fit parameters are 
given. We choose Ho = 50kms~^ Mpc~^ . 
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Table 4. Evolution characteristic times: look-back time to 
which all luminosities were higher by a factor e, given the 
PWLE parameter in Table 2. Times are in units of H^^ , the 
mean is 0.27 and the dispersion is 0.04. 



Table 6. PWLE characteristic times (Table 4) in units of the 
age of the universe in the corresponding cosmology. Mean: 0.31; 
dispersion: 0.02. These ratios are far more independent of cos- 
mology than the times in Table 4. 
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Table 5. Same as 


Table 4, but for LEXP. Mean characteristic 
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time: 0.13; dispersion: 0.05. 
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Table 7. Same as Table 6, but for LEXP. Mean; 0.15; disper- 
sion: 0.02. As for PWLE, the characteristic time depends far 
less on the age of the universe than on the Hubble time. 
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